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Abstract 



We study the bound-free electron-positron pair production in relativistic 
heavy ion collisions to different bound states in a full Plane Wave Born Ap- 
proximation (PWBA) calculation. Exact Dirac wave functions are used for 
both the bound electron and the free positron in the final states. Results for 
the Relativistic Heavy Ion Collider RHIC as well as the forthcoming Large 
Hadron Collider (LHC) are given. This process is one of the dominant beam 
loss processes and can become critical for the operation with heavy ions. A 
simple parameterization is given as well. We compare our results with calcu- 
lations of other groups. 

I. INTRODUCTION 

Bound free pair production is one of the new types of processes that occur in relativistic 
collisions of atoms and ions. It is the production of an electron-positron pair with the 
electron not produced as a free state but as a bound state of one of the ions: 

Z P + Z T -> Z P + (Z T + e-) + e+. (1) 

This process changes the charge state of the ion. Due to the change in the charge-to-mass 
ratio, such an ion will be lost from the circulating beam and the luminosity will be seriously 
affected. Together with the electromagnetic dissociation this process is the dominant beam 
loss process at RHIC and the Large Hadron Collider (LHC) at CERN when using heavy ions 
|T|-[3[|. In addition to the beam loss itself, it was recently discussed in [|| that the capture 
process can lead to localized beampipe heating. This could cause magnet quenches if the 
local cooling is inadequate. Therefore it is very important to know this cross section with 
high accuracy and reliability. 

We want to mention also that the corresponding process using antiprotons as the "target" 
has recently found an application in the production and detection of relativistic antihydro- 
gen, see [j5]|| . In a collision of the antiprotons with a target again electron positron pairs are 
produced. The positron can be produced in a bound state of the antiproton forming an anti- 
hydrogen atom. Using the equivalent photon method, this cross section was calculated in [|7| . 
In the equivalent photon approximation, the cross section factorizes into a photo-induced 
cross section and an equivalent photon spectrum. This spectrum depends on a cut-off pa- 
rameter which should be chosen appropriately. Using the PWBA expression for the cross 
section one can avoid this ambiguity and see also how this parameter should be chosen. This 
was done in || and ||. In these papers, differential as well as total cross sections are given 
as a function of the collision energy. Whereas in the latter reference approximate lepton 
wave functions appropriate for low values of Zt, the charge of the antiproton, have been 
used, exact Dirac wave functions, which are valid also for higher values of Zt, have been 
used in ||. 

It is the purpose of this paper to calculate the cross-sections for the energy region of 
the colliders RHIC and LHC exactly within the plane wave Born approximation (PWBA) 
or (which is equivalent) the semiclassical straight line approximation (SCA), including also 
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higher shells. Previously, such kinds of calculations were done using the Weizsacker- Williams 
method [10JL1]. Higher order effects (in the interaction with the projectile) have been 
considered by a number of groups, see, for example, [p^-p^|. Recently exact calculations 
have been done in the high energy limit by Baltz ||15|| . He finds that that the contributions 



from higher orders are rather small (of the 1% level) and even tend to decrease the cross 
section. This is in contrast to results found at much smaller energies with Lorentz factors 
in the target rest frame 7t < 3 fT^JT^l > but see a ^ so ■ 

In Sec. |IJ we extend the formalism presented in [|J to the relativistic heavy ion case, 
including now also formulae for the higher atomic states. Numerical results are then pre- 
sented in Sec. |T1| and we discuss the dependence on the beam energy and on the principle 
quantum numbers n and k. We review these results with those existing in the literature in 
Sec. [TV|. Finally our conclusions are given in Sec. [V]. 



II. PWBA OR SCA THEORY OF BOUND-FREE PAIR PRODUCTION 



The total cross section for bound-free pair production (per electron-state) in Lorentz 
gauge is given by p0| , 
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where (3t = - is the velocity of the projectile in the target rest frame, 7 T = (1 — /3f>) is the 
Lorentz factor in the target rest frame, and F(k 2 ) is the form factor of the charge distribution 
of the projectile. For electrons and positrons we can set this form factor F(k 2 ) = 1. (For the 
pair production of muons and tauons this formfactor will become important.) We denote 
the fine structure constant by a em to distinguish it from the usual Dirac matrices a. The 
spatial momentum transfer from projectile to target is hk . The absolute value of the wave 
vector k is denoted by k, and its z-component is related to the energy of the photon as 
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The total energy of the bound electron is Ef and the one of the electron in the continuum 
is Ei. (Please note that Ei is negative.) ipi(r) and ipf(r) are the Dirac-Coulomb wave 
functions describing the initial negative continuum state and final bound state. The charge 
numbers of the projectile and target are denoted by Zp and Zt- 
The same expression Eq. (H) in the Coulomb gauge reads 
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Eq. (Q) can be obtained from Eq. @ by making use of the current conservation. Both 
equations are therefore identical only when exact eigenfunction of the Dirac Hamiltonian 
are used, as already stressed in pj| , [2T|| . 

The 4-component Dirac spinor in a spherically symmetric field is given by 



ifK.(r)x™{r) 



The angular dependence is expressed by the spin-angular functions 
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where \ T is the Pauli spinor and 
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The radial Dirac equation is 
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The radial functions in the Coulomb field V(r) = — £ (with £ = a em Zy) for the bound states 
are given by 
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The corresponding energy eigenvalues are 



E n j = mc 



1 + 



n-j-± + 7. 
where n is the principle quantum number defined by 

n = n r + \k\. 

The quantity (3 corresponding to the energy eigenvalue is 
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7 is given by 
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The adequate normalization for bound states is 



d 3 r ip'(r)ip(r) = I dr r A 
o 



and the normalization constant is found to be 
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which are normalized according to 
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In ( |TT| ) k', i], ip and Nf are given by 

(E 



k' 



E 2 - (mc 2 ) 2 

he 



tick' 



7 — if] 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



2 7 e^|r(7 + l + ^)| 

r(2 7 + i) 



(20) 



Please note that E > mc 2 for positive energy continuum wave functions and E < —mc 2 for 
negative energy states. 

Starting from the expression in the Coulomb gauge (|j) the angular integration can be 
performed analytically |22] to get: 
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a (bfpp) = lg 2 
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In this form the increase of the cross section with In 7^ can be clearly seen due to the second 
term (proportional to Tj_) in the integral. T\ and T± are given by 
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and the radial integrals are 
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These rapidly oscillating integrals are very difficult to evaluate. The integrals T} and Tj_ 
are independent of the Lorentz factor 7 and can be evaluated before integrating over k in 
Eq. (|2l"l). But for high values of 7 this integration extends from very low to very high values 
of k, as well as, Ei. We were able to do the integration for these extreme limits for RHIC 
and LHC For this we have used recursion relations from 



III. NUMERICAL RESULTS 

We give our results for symmetrical collisions where the charges of both ions are equal 
Zp = Zt = Z . The cross sections refer to capture to one of the ions only. Since the cross 
section for capture scales with Zp, the cross section for asymmetrical collisions is obtained 
by scaling the cross sections with Zp/Z 2 . 



6 



In Table | we give the cross sections for the different ions and for different bound states. 
This table was already presented in [22,24]. We have used the conditions as they are relevant 
for Au-Au at RHIC and Pb-Pb at LHC In these colliders each of the ion is assumed to have a 
Lorentz-factor of 7 C = Z/A'-fp. For RHIC we take 7 P = 250, for LHC 7 P = 7500, respectively. 
For different values of 7 C the cross section can be obtained by using Eq. (p8|). The cross 
section are larger than the ones quoted in the 'ALICE technical report" [|3|. The Lorentz 
factor 7t in the rest frame of one of the ions, which is the relevant one for the bound-free 
pair production, is then given by 



7t = 2 76 2 - 1. 

For RHIC (LHC) this corresponds to 7 T = 2 x 10 4 (7 T = 1.8 x 10 7 



it was found that the cross sections can be parameterized as [13 



a = A ln(7 c ) + B 



(27) 

For large values of 7 C 
(28) 



We give also the values of the constants A and B in the table. 

In Figure [1] we show the dependence of the different cross sections on the ion-charge Z. 
As the cross section for the s states is known to be proportional to Z 7 for low values of Z 
2~5| , we divide through this factor. The most striking result is that the cross section for 
the capture into the 2pi/ 2 state increases rapidly with Z and becomes even larger than the 
one for the 3s state for heavy ions. The cross section for the 2p 3 / 2 shows a clearly different 
behavior. We attribute this difference to the "small component" of the Dirac wave function. 
This wave function has an s-wave character for the pi/2 state and a <i-wave character for the 
P3/2 state, respectively. With the increase of Z the "small component" gets more important 
and leads to the different behavior of the two cross sections. This is also in accord with a 
qualitative discussion of relativistic effects of the Dirac wave function for bound states as a 
function of Z |26| . 

For the cross section to the different s states as a function of the principle quantum 
number n, a scaling law is found in the nonrelativistic limit, 



n 3 



(29) 



that also appears, e.g., in the photo-electric effect |26[. This scaling law arises because the 
value of the wave function at the origin 



(30) 



enters in the expression for the cross section. Our calculations confirm this 1/n 3 scaling to 
a high accuracy once more and for all values of Z. With this scaling law we can sum the 



contribution from capture processes into all s states [25 



a ns = C(3) vis ~ 1-202 a ls 



(31) 



n=l 



where ( denotes the Riemann (^-function. For low atomic numbers Z the contributions from 
s-orbits are the most important ones, whereas the cross section to the 2pi/2 contributes of 
the order of 6.5% in the heaviest cases. 
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IV. COMPARISON WITH OTHER CALCULATIONS 



In this section we want to give a comparison of the results that have been obtained with 
different approaches. The results for the cross section are summarized in Table [[]]■ A large 
number of calculations exist including the electromagnetic interaction (with the projectile) 
to all orders. Most of them have been done for beam energies much smaller than those 
available at the relativistic heavy ion colliders. In order to get results for energies at RHIC 
and especially at LHC energies extrapolations were used, see below. 



A calculation within the PWBA approximation was done by Becker et al. |27].|28].[T4 
Calculations up to values of 7^ = 1000 have been done and an interpolation formulae of the 
form 

a « Z\Z\a In ( Tt / 7 o) (32) 



(using our notation) and tabulated values for a and 70 can be found in ||27| . The parameters 



a and 70 depend slightly on Zt- Using the value of a for Z = 80 for both Au and Pb, we 
get values of 93 barn (Au-Au at RHIC) and 226 barn (Pb-Pb at LHC), respectively. 



In [£5| a calculation is done making use of the Sommerfeld-Maue wave functions as an 
approximation to the exact Coulomb-Dirac continuum wave functions and using also an 
approximation for the bound state wave function, which is valid in the limit Za <C 1. 
They found their result in SCA to be equal to the one derived within the equivalent photon 
approximation. Using their formula we obtain results which are substantially lower (by a 
factor of two) than our results here. As their approximation is strictly valid only in the limit 
Za <C 1, such a discrepancy is perhaps not surprising. 

Baltz et al. treated the problem in a series of papers |12,13,|T8,29|| . For large impact 



parameter they also make use of a perturbative treatment. In addition they simplify the 
interaction potential by taking only lowest order terms in 1/ 7 t and also expanding in p/b, 
where p is the electron coordinate and b the impact parameter. They proposed the param- 
eterization of the cross section for large values of jt of the form 

a = Aln(j T ) +B, (33) 

with the interpretation of A to be given by the perturbative part only and the influence 
of higher order terms at small impact parameter to be present in B only (Of course this 
parameterization is identical to the one of Becker et al. above). Tabulated values for A 
and B are given in |13|| . From them (and choosing for b m i n = 1A C ) we get cross section of 
83 barn (Au-Au at RHIC), 161 barn (Au-Au at LHC), and 436 barn (U-U at LHC). They 
also give the cross section to be 89 barn (Au-Au at RHIC) in the text where contributions 
from nonperturbative processes at small impact parameter have been included. No results 
for lead-lead collisions are given, but assuming a Z 7 scaling we get from the extrapolation of 
the Au results a value of 206 barn (Pb-Pb at LHC), from the extrapolation of the U results 
195 barn (Pb-Pb at LHC) using the jt dependence of Eq. (|33|). 

In [|30| a calculation is done using the cross section of the free pair production cross 
section and folding it with the momentum distribution of the bound state. The electron 
therefore is described by a plane wave. For Au-Au collisions at RHIC a cross section of 
72 barn was found by them, but see also [pTJ. 
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Two calculations making use of the equivalent photon approximations have been per- 
formed also. In both exact solutions of the Dirac equations are used. In the equivalent 
photon method a cutoff parameter has to be introduced, which leads to some uncertainty 



in the results. In |L0| capture to the ls-state was calculated. The results were also com- 
pared with the cross section for bound-free pair production induced by real photons at low 
photon energies by |32 |. For the heavy ion case values of 90 barn (Au-Au at RHIC) and 
222 barn (Pb-Pb at LHC, 7 C = 3400) are given, respectively. Also tabulated values of A 
and B are presented. Their results are in good agreement with the present values. The 
cutoff parameter was taken to be the Compton wavelength of the electron. In [IT] similar 
calculations are performed including also capture cross section calculations to the L-shell. 
Agreement with the results of fH| was found for the .fT-shell. A critical discussion of the 
use of the cutoff-parameter present in the equivalent photon spectrum was done. The ex- 
act value of this parameter is not given by the Weizsacker- Williams theory and therefore 
introduces some uncertainties in the results. For Au-U collisions at RHIC with the electron 
being captured by the U-atom they find a value of 165 barn for b m i n = 2A C and 182 barn for 
bmin = A c with A c the Compton wavelength of the electron. Comparing these results with 
the 195 barn we get, an cutoff parameter of b m i n = A c seems to be favored. Again assuming 
a scaling of the result with Z\ we get for Au-Au collisions at RHIC 77 barn {b min = 2A C ), 
85 barn (b m i n = A c ), respectively. Unfortunately no results for LHC energies are given. 

Experiments were first done at Bevalac at 1 GeV/A [j33] -|35"ll . A cross section of 
2.19(0.25) barn for Au-U collisions and capture on the U projectile was found. Experi- 
ments were also done at the AGS at 10.8 GeV/A corresponding to 7t = 12.6, 7 C = 2.6 



36| , |3"7j| . A cross section of 8.8 barn was found for Au-Au collisions (Our calculation gives 



11.86 barnQ The experimental results were found to be in agreement with the theoretical 
results of ©[50]). 

A measurement has recently been done at the CERN SPS with a 158 GeV/A Pb beam 
(jt ~ 168, corresponding to 7 C = 9.2) EH], see also pO]. The capture to the Pb ion was 



studied for several different targets. For a gold target a value of 44.3 barn is given (our 
result is 45 barn). Assuming a scaling of the cross section as given above Eq. ( ^8|) and 



assuming B = —24 barn from fI8|, they give extrapolated values of 94 barn (Au at RHIC) 



and 204 barn (Pb at LHC). These values again agree quite well with our results. 



V. CONCLUSIONS AND OUTLOOK 

Apart from neglecting higher order effects in the projectile charge, we have given a 
full calculation of the bound-free capture cross-section. Since such higher order effects 



were shown to be small ||15|| , our predictions should be very reliable. Overall our results 
are in agreement with a number of other calculations, using different approaches. These 
calculations are also very important for the questions of luminosity loss and localized beam 



Please note that while we compare cross section for the ls-state only in the comparison with 
other calculations, comparing with experimental results we use the total cross section into all bound 
states. 
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pipe heating [f|, especially at the LHC. 

Muon and r-lepton bound-free pair-production can be calculated similarly. Since the 
bound-free pair-production cross section scales with l/mf epton |25j these cross-sections will 



be much smaller than the ones for e + -e~ pair-production. Especially for heavy nuclei, it will 
be necessary to use finite size lepton wave functions. In addition to the l/mf epton scaling 
there will be also a severe reduction of the cross section due to the formfactor effects. Maybe 
such exotic atoms could be interesting for physics JIIJ? We mention that the lifetime of the 
free tau-lepton is given by ct = 87.2//. For a Lorentz factor of about 100 this is still a very 
short distance of about 0.1mm for the decay length of such an atom. On the other hand 
a muonic atom lives long enough and could be extracted from the beam (similar to the 
antihydrogen experiments |5|||). 

In a fixed target experiment, the effect of screening of the atomic electrons can be incor- 
porated as well by using a screened Coulomb potential for both the wave function, as well 
as, the interaction with the projectile (F ^ 1 in Eq. (Q)). This could be relevant for the 



CERN fixed target experiments This effect has been studied in and was found 

to be on the percent level. 

RHIC is running right now and we expect that the present numbers will soon be tested. 
After the preparation of this manuscript, an article by Bertulani and Dolci appeared E2 



This is a continuation of H to the heavy ion case. In it they give an explanation why the 



approximation used in [25] does not work well for large values of Z. 
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FIG. 1. The cross section for the capture of the electron into different bound states is given as 
a function of the charge of the ion charge Z. A Lorentz factor in the cm. system of 7 C = 3400, 
corresponding to an equivalent Lorentz factor of jt = 2.3 x 10 7 in the rest frame of the ion is used. 
For s-states there is an approximate Z 7 scaling of the cross section. 
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TABLES 



1 1 , 1 

bound state 


cj(RHIC) 


\1 1 

[b] 


<x(LHC)[&] 


A[b] 






B[b] 






1 TT 1 TT 


7c = 


250 




7c = 


7500 














Is 


2.62 


• 10" 


1 1 


4.25 


io- n 


5.36 


10" 


1 9 

- 1Z 


-3.40 


10" 


1 9 
-1Z 


2s 


3.28 


• 10" 


1 9 

-1Z 


5.31 


KT 12 


6.70 


io- 


1 '1 

- 11 


-4.23 


10" 


1 ,; t 

- LO 


2p(l/2) 


3.75 


• 10" 


17 


6.10 


io- 17 


7.73 


10" 


-18 


-5.20 


10" 


-18 


2p(3/2) 


1.47 


• 10" 


17 


2.41 


io- 17 


3.10 


10" 


-18 


-2.42 


10" 


-18 


3s 


9.70 


• io- 


13 


1.57 


io- 12 


1.98 


io- 


-13 


— 1.26 


10 


-13 


/u Ca- u Ca 


7c = 


125 




7c = 


3750 














Is 


1.61 


• io- 


_2 


2.92 


10~ 2 


3.84 


10" 


_3 


-2.48 


10" 


_3 


2s 


2.00 


• 10" 


Q 
O 


3.62 


10~ 3 


4.78 


10" 


_ 1 


-3.07 


10" 


„4 


2p(l/2) 


1.39 


• 10" 


_5 


2.52 


10~ 5 


3.35 


io- 


_5 


-2.33 


10" 


_5 


2p(3/2) 


3.63 


• io- 


-5 


6.70 


10~ 6 


9.02 


10" 


_7 


-7.27 


10" 


-7 


3s 


5.90 


• 10" 


-4 


1.07 


10" 3 


1.41 


10" 


-4 


—9.10 


10" 


-5 


A7 a 47 a 

4 'Ag- 4 'Ag 


7c = 


109 




7c = 


3264 














Is 


3.51 






6.46 




8.68 


10" 


_1 


-5.63 


10" 


_1 


2 s 


4.33 


• io- 


i 

- 1 


7.98 


io- 1 


1.07 


10" 


1 

- 1 


-6.94 


10" 


o 
-A 


2p(l/2) 


2.81 


• 10" 


-2 


5.21 


IO" 2 


7.05 


10" 


_3 


-5.02 


10" 


_3 


2p(3/2) 


3.80 


• 10" 


.3 


7.16 


10" 3 


9.87 


10" 


_4 


-8.31 


10" 


_4 


3s 


1.26 


• io- 


-1 


2.34 


io- 1 


3.13 


10" 


-2 


—2.02 


10" 


-2 


7Q a 7Q a 


7c = 


100 




7c = 


3008 














Is 


94.9 






176 




23.8 






-14.7 






2 s 


12.1 






22.4 




3.04 






-1.87 






2p(l/2) 


3.62 






6.77 




9.27 


10" 


-1 


-6.56 


10" 


-1 


2p(3/2) 


2.10 


• io- 


x 


4.01 


io- 1 


5.62 


10" 


-2 


-4.93 


10" 


-2 


3s 


3.46 






6.40 




8.67 


10" 


-1 


—5.34 


10" 


-1 


8 "Tb- 2 Pb 


7c = 


99 




7c = 


2957 














Is 


121 






225 




30.4 






-18.7 






2 s 


15.5 






28.8 




3.91 






-2.39 






2p(l/2) 


5.21 






9.76 




1.34 






—9.46 


10" 


-1 


2p(3/2) 


2.78 


• 10" 


-1 


5.33 


1 n— 1 

10 


7.50 


10" 


-2 


-6.61 


10" 


-2 


3s 


4.42 






8.20 




1.11 






-6.79 


10" 


-1 


92 U- 92 U 


7c = 


97 




7c = 


2900 














Is 


263 






488 




66.0 






-39.0 






2s 


34.4 






63.7 




8.63 






-5.10 






2p(l/2) 


16.7 






31.3 




4.30 






-3.00 






2p(3/2) 


6.77 


• 10- 


1 


1.30 




1.83 


10" 


-1 


-1.63 


10" 


-1 


3s 


9.67 






17.9 




2.43 






-1.44 







TABLE I. Cross section for the bound-free pair production of one ion only for different bound 
states are given for RHIC and LHC conditions for different ion-ion collisions. Also given are the 
parameters A and B to be used in Eq. ( p8[ ) for the dependence on the Lorentz factor j c . 
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a(Au-Au,RHIC,7 c = 100) [barn] 
94.9 



CT(Pb-Pb,LHC,7 c = 2957) [barn] 
225 



Ref. 
This work 



93 


226 




E 


37 


86 


2| 


89 


206 (from U) 


T|l3,18] 




195 (from Au) 






72 




a 




90 


222 ( 7c = 3400) 


m 





77 (6, 
85 



2A C ) 

'min — ^c) 

94 



204 barn 



[11] 



TABLE II. Cross section for bound-free pair production of one ion only into the ls-state 
calculated by different groups are compared. The last entry uses experimental results of bound-free 
pair production into any bound state at CERN SPS with 158 GeV/A energies. Given are the results 
for Au-Au collisions at RHIC and Pb-Pb collisions at LHC. For details of the different calculations 
see the text. 
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